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The reaction mechanism of the central collisions and peripheral collisions for 
112, 124 Sn + H2,124 Sn at = §0MeV is investigated within the framework of the 

Improved Quantum Molecular Dynamics model. The results show that multifragmenta- 
tion process is an important mechanism at this energy region, and the influence of the 
cluster emission on the double n/p ratios and the isospin transport ratio are important. 
Furthermore, three observables, double n/p ratios, isospin diffusion and the rapidity 
distribution of the ratio R 7 for 112 - 124 Sn +112,124 Sn at E/A=50MeV are analyzed 
with the Improved Quantum Molecular Dynamics model. The results show that these 
three observables are sensitive to the density dependence of the symmetry energy. By 
comparing teh calculation results to the data, the consistent constraint on the density 
dependence of the symmetry energy from these three observables is obtained. 

1. Introduction 

The Equation of State (EOS) of asymmetric nuclear matter can be written approx- 
imately as E(p, 8) =E(p,S = 0) + E sym (p)6 2 + 0(5 4 ) , with 6 = (p n -Pp)/(pn + Pp)- 
It closely relates to different areas of nuclear physics, such as the nuclear structure 
of finite nuclei, dy namics proc ess of neutron rich heavy ion collisions, physics of 
neutron star, et.alP^ ^^®^ . For the symmetry nuclear matter, measurements 
of isoscalar collective vibrations, collective flow and kaon production in energetic 
nucleus-nucleus collisio ns have constrained the EOS from normal density to five 
times saturation density^ ^ l lO l However, large uncertainties exist on the theoret- 
ical prediction on the density dependence of the EOS for neutron matter ^H. So, 
determining the density dependence of the symmetry energy, E sym (p), becomes 
one of the main goals in nuclear physics at present and in the near future and has 
stimulated many theoretical and experimental studies. Heavy ion collisions with 
neutron-rich nuclei provide a unique opportunity to obtain the information of the 
density dependence of the symmetry energy in the laboratories because large extent 
of density can be formed d uring the HICs. Many usef ul observabl es from HICs, such 
as isoscaling of the isotope EEHH^ j S0S p m diff us i on | 15 | 16 | 17 118^ neutron to proton 
yield ratios and its flow at intermediate heavy ion collisions have been proposed 
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and studied both theoretically and experimentally in order to constrain the density 
dependence of symmetry energy at subsaturation density. And n~ /n + ratios and 
its flow in high energy HICs was proposed to c onstrain the density dep endence of 
the symmetry energy at supersaturation density^ 2 ^ 21 ! 22 ! 2 ^ 24 ! 2 ^^ ! 

Since the complex cluster production are observed in experiment for interme- 
diate energy HICs, it is necessary to investigate the effects of clusters emission on 
these observables, neutron to proton ratios and isospin diffusion. As a microscopic 
dynamical n-body transport theory, formation of fragments is included automati- 
cally in the QMD type model. So, the ImQMDO J EP7 | 28 | 34 | 

model is suitable to 

study the impacts of the cluster formation on the isospin sensitive observables, 
such as n/p ratios for emitted nucleons and isospin diffusion for central collision 
and peripheral collisions at E/A — 5QMeV for Sn + Sn. Fig.l shows the contour 
plots of multiplicity of fragments with charge Z as a function of its scaled rapidity 
v/Vblam at E/A = 50MeV calculated with soft symmetry energy (ji = 0.5) for 
different impact parameters, b = 2, 6, 8/m, respectively. In general, one can observe 
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Fig. 1. The contour plots of multiplicity of fragments with charge Z as a function of its scaled 
rapidity y/y^ m for 112 Sn + 112 Sn at E/A = 50MeV for b = 2, 6, 8/m with 7, = 0.5 cases. 



two peaks on the rapidity distribution for the charge Z of fragments from cen- 
tral collisions to peripheral collisions, and the intermediate mass fragments can be 
formed from middle rapidity to projectile/target rapidity for central collisions and 
peripheral collisions. In detail, the maximum charge of fragments at mid-rapidity 
decrease from Z ~ 15 to Z ~ 8 with the impact parameter b increasing, and same 
for its yield. Furthermore, the heaviest fragments have lost about 35% of their ini- 
tial velocity for central collisions and about 25% (10%) of their initial velocity for 
b = 6/771 (b = 8/m). The quantities of velocity loss for heaviest fragments depend 
on the decelerating effects from the effective n-n interactions and nucleon-nucleon 
collision in the participant region, and the multiplicity of fragments formed below 
the normal density in HICs. Through above discussion, one can see that the clus- 
ters formation plays an important role at intermediate energy HICs from central 
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to peripheral collisions and the effects of cluster emission on the isospin sensitive 
observables could not be ignored. 

In order to study the density dependence of symmetry energy and the influence 
of cluster emission on the isospin sensitive observables, the density dependence of 
symmetry energy in nuclear matter: 

TTi / \ 1 ft2 2/3/37T 2 P ,2/3, , / n \ 

6 2m 2 po 2 po 

is introduced. C s is the symmetry potential strength parameters, and ji give the 
density dependence of symmetry energy. The different density dependence of sym- 
metry energy can be realized by varying ji in the ImQMD model. 



2. The n/p ratios, isospin diffusion and the cluster emission effects 

The neutron to proton ratio R n / V = ~jjgr^ / ^ °^ P r e-equilibrium emitted 
neutron over proton spectra was considered as a sensitive observable to the density 
dependence of symmetry energy^} because it has a straightforward link to the 
symmetry energy. In order to reduce the sensitivity to uncertainties in the neutron 
detection efficiencies and sensitivity to relative uncertainties in energy calibrations 
of neutrons and protons, the double ratio 

DR(n/p) = R n/ JA)/R n/ JB) = dM n (A)/dE cm dM n (B)/dE cm 

" M " n/pK ' dM p {A)/dE c J dM p {B)/dE cm K) 

had bee n mea sured by Famiano and compared with the transport model 
predictio n 1 2 ^ 21 !. 

We have performed calculations of collisions at an impact parameter of b = 2/m 
at an incident energy of 50AIeV per nucleon for two systems: A = 124 Sn -f 124 Sn 
and B = 112 Sn + 112 Svfil with ImQMD05 to study the DR(n/p) ratio for 
emitted nucleons. The shaded regions in the left panel of Fig. 2 (a) show the 
range, determined by uncertainties in the simulations, of predicted double ratios 
DR(n/p) — R n / p (A) / R n / p (B) of the nucleons emitted between 70° and 110° in the 
center of mass frame as a function of the center of mass nucleon energy, for 7,; = 0.5 
and 2. The double ratios DR(n/p) are higher for the EOS with the weaker symme- 
try energy density dependence 7, = 0.5 than that for 7,; = 2.0 because the nucleons 
mainly emit from the lower density region at intermediate energy HICs. Compare 
to the data on DR(n/p) for emitted nucleons (solid stars), the general trend of data 
DR(n/p) are qualitatively reproduced and the data seem to be closer to the calcula- 
tion employing the EOS with 7; = 0.5. The fig. 2 (b) show the coalescence-invariant 
double ratio. The coalescence-invariant double ratios are constructed by including 
all neutrons and protons emitted at a given velocity, regardless of whether they 
are emitted free or within a cluster. The data are shown as open stars and the 
calculation results are shown as shaded region in the right panel of figure 2(b). 
Here, the measurements and simulations results illustrate that the fragments with 



November 9, 2009 9:4 WSPC/INSTRUCTION FILE iwnd09zhyx 



4 Yingxun Zhang 

Z > 2 mainly contribute to the low energy spectra and do not affect the high-energy 
DR(n/p) data very much. 




Fig. 2. (a) DR(n/p) ratios for emitted free nucleons and (b) coalescent-invariant DR(n/p) from the 
ImQMD simulations are plotted as shadow region. Data (star) from NSCL 19 .(c) DR(n/p) ratios 
for d iffer ent ji for central collisions, Data (star) from NSCL. (d) x 2 analysis on the DR(n/p) 
ratioJSSl. 

In order to constrain the range of ji from the DR(n/p) data that had been 
published, a series calculations for two systems, A = 124 Sn + 124 Sn and B = 112 
Sn + 112 Sn, have been performed by varying ji — 0.35,0.5,0.75,1.0 and 2.0(linc 
with symbols j^Sl. As shown in fig. 2(c), the computation uncertainties are statis- 
tical. It is known that emitted nucleons mainly from the subnormal densities at 
this energies. So, the n/p ratios of emitted nucleons are associated with the values 
of symmetry energy at subnormal density. Therefor, the DR(n/p) ratio should in- 
crease with decreasing 7.;. However, in the limit of very small 7* -C 0.35, the system 
completely disintegrates and the DR{n/p) ratio decrease and approaches the limit 
of reaction system, [N / Z)^/ {N / Z)n2 — 1.2. As a consequence of these two com- 
peting effects, the double ratio values peak around the 7$ = 0.7. Despite the large 
experiment uncertainties for higher energy data, those comparisons definitely rule 
out very soft (7$ = 0.35) and very stiff (7; = 2.0) density dependence of symmetry 
energy. Fig. 2(d) shows the dependence on 7, of the total \ 2 computed from the 
difference between predicted and measured double ratios. Within a 2er uncertainty, 
parameters of 7^ fall in the range of 0.4 < 7; < 1.05 for the C s — 35.2MeV. 

It has been proved that the isospin d iffusion ability depends on the strength of 
the symmetry energy in HICsISIIEES] i n order to quantify the isospin diffusion 
degree in heavy ion collisions, the isospin transport ratios Ri 

Ri = (2X — Xa+a — Xb+b)/(Xa+a — Xb+b) (3) 

has been introduced^. The subscript A and B represent the neutron rich and 
neutron-poor nuclei, and the A = 124 Sn, B - 112 Sn in this work. The value of Ri 
is obtained through 3 reaction systems at least, A + A, B + B and A + B(orB + A) . 
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Where X is the isospin tracer from the isospin asymmetry nuclear reaction sys- 
tem A + B(orB + A). The non-isospin diffusion effects are minimized by scaling 
the isospin observables with the same observables from the symmetric collisions 
of the neutron-rich A + A and neutron-deficient B + B systems using the isospin 
transport ratio defined as above. Adopting the Eq.3 definition on Ri, one can ex- 
pects Ri = ±1 in absence of isospin diffusion. In the opposite, Ri « if isospin 
equilibrium is achieved. Eq.3 also dictates that d ifferent observable, X, can give 
the same results if they are linearly relatecP^^A In experiment, the X is taken 
as the isoscaling parameter X a = a and th e yield ratios of A = 7 mirror nuclei, 
X 7 = ln(Y( 7 Li)/Y( 7 Be)) for peripheral HIC J 15 | 17 1 In order to analyze the isospin 
diffusion data with transport model, the observable X = 6, the isospin asymme- 
try constructed from the fragments and free nucleons at the relevant rapidities is 
adopted in ImQMD. It has been confirmed theoretically and experimentally that 
there is a linear relationship between X exp = X7, X„ and Xth. = 5, and the rela- 
tionship R(X = a) = R 7 {X = X 7 ) = R{6) hold J 13 1 32 ' 33 - 13 ^ 1 . 

The left panel in Fig. 3 show the ImQMD predictions on the R(X — 8) (lines) 
near the projectile ra pidi ty as a function of impact parameters b for different ji = 
0.35, 0.5, 0.75, 1 and 2^. Faster equilibration occurs for smaller 7* values which 
correspond to a larger symmetry energy at subnormal densities. Thus, we see a 
monotonic de crea se of the absolute values of Ri(S) with decreasing 7$. Experimental 
data on Ri(a}^, is plotted as shaded regions in the left panel of Fig. 3. Performing 
the x 2 analysis between the ImQMD predictions and the experimental data on 
Ri(a), the range of the symmetry potential parameter 0.45 < ji < 1.0 are obtained 
within 2cr. Furthermore, the ImQMD predictions on the R(X = 5) (lines) as a 
function of rapidity are calculated and shown in the middle panel of fig. 3 as lines, 
the star symbols represent measured values of Rj obtained from the yield ratios of 
7 Li and 7 Be at b = 6/rJ3. It can be seen that the calculation of the shapes and 
magnitude of the rapidity dependence of the isospin transport ratios reproduces the 
trends accurately. The corresponding x 2 analysis on Rj with calculations at b = 6 
and 7/m favors the region 0.45 < 7, < 0.95 within 2a uncertainties. 

In order to see the cluster effects on the isospin transport ratio Ri, a tracer 
defined by the isospin asymmetry of the heaviest fragments with Z > 20 in projectile 
region is adopted in the transport model simulation for peripheral HICs (b = 6/m). 
The right panel in fig. 3 show the Ri calculated with X = S(Z max > 20) (open 
diamonds) and open circles represent the results from the X — S near the projectile 
rapidity as previous discussed. Independent of the isospin tracer we adopted, Ri 
obtained with soft symmetry case (7$ = 0.5) is smaller than those obtained with 
stiff symmetry potential case (7^ = 2.0). Another side, the values of Ri obtained 
with the new tracer X = S(Z max > 20) are larger than that with X = 5 constructed 
from emitted nucleons and fragments near the projectile rapidity. It is caused by 
the cluster emission effects and the dynamical properties of the isospin diffusions. 
In HICs process, the neutrons (protons) of the projectile or target mainly transfer 
to the lower density neck region and then break up into IMFs and nucleons, few 
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Fig. 3. Left panel: The calculated results on isospin transport ratios from the ImQMD model(lines) 
as a function of impact parameters for different values of fi and the experiment data Ri(a) 
(shaded regions). Middle panel: ImQMD calculations on the isospin transport ratios as a function 
of rapidity fo r di fferent 74 at b = 6/m and the experiment data on R7 as a function of rapidity 
(star symbols jSSl Right panel: Isospin transport ratios obtained from different isospin tracer, X = 
S(N, Frags)) (open circles) and X = 8(Z max > 20) (open diamonds) in projectile region 

nucleons from target residues diffuse to the projectile residues or vice versa due to 
the short diffusion time scale at E/A — bQMeV. It leads the larger values of Ri 
for the new tracer X — 6(Z max > 20). So, study the Ri with these different isospin 
tracer will help us to understand the dynamic effects on isospin diffusion in HICs 
and give a new constraint on the density dependence of symmetry energy. 

3. Summary and outlook 

In summary, the reaction mechanism of the central and peripheral collisions are 
investigated for 112 ' 124 SVi + 112 . 124 Sn at E/A = 50MeV for different impact pa- 
rameters. The results show that the reactions systems do not reach the global 
thermal equilibrium even for Sn + Sn central collisions at E/A = 50MeV and the 
effects of cluster emission are important for both central and peripheral collisions. 
Three observables, double n/p ratios, isospin diffusion and the rapidity distribu- 
tion the ratio R7 are analyzed with the Improved Quantum Molecular Dynamics 
model. The results show that these three observables are sensitive to the density 
dependence of symmetry energy. By comparing calculation results to the data, the 
consistent constraints on the density dependence of symmetry energy can be ob- 
tained. Furthermore, the analysis of the coalescence-invariant DR(n/p) ratio and 
the isospin transport ratios Ri(X — 6(Z max > 20)) show that the influence of the 
cluster emission on the double n/p ratios and the isospin transport ratio are impor- 
tant. However, the large uncertainties on the data implies that the improvements in 
the precision of data should be done in the near future, and accurate constraints on 
the density dependence of symmetry energy should be complemented from many 
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aspects, such as giant dipole resonances, pygmy dipole resonance, and mass data 
and some new observables from HICs and the observed properties of neutron star. 
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